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基于 光合 有 效 辐射 瞬时 值 估算 日 均值 的 方法 ' 
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摘 要 : 光合 有 效 辐 射 是 太阳 辐射 的 一 部 分 ,是 植被 进行 光合 作用 和 陆地 
生产 力 具 有 重要 影响 。 本 文 基于 SURFRAD 多 年 观测 数据 ,建立 了 一 个 线性 回归 模型 ,由 上 午 、 下 午 瞬 时 光合 有 效 


= 


生态 系统 碳 循环 的 核心 因素 ,对 估算 植被 


辐射 值 估算 日 均值 。 该 模型 模拟 效果 较 好 ,单个 观测 站 的 均 方 误差 低 于 9 W m? ,判定 系数 不 低 于 0.96, 并 适用 
于 不 同 经 纬度 、 气 候 条 件 和 海拔 高 度 下 的 区 域 ,总 的 均 方 误差 为 8.863 1 Wem”, 判定 系数 R 为 0.977 ,表明 该 模 


型 和 方法 有 较 好 的 可 行 性 和 可 靠 性 。 
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光合 有 效 辐射 (PAR ，photosynthetically active 
radiation ) 是 太阳 辐射 在 电磁 波谱 的 可 见 光 (400 ~ 
700 nm) 部 分 ,是 植被 进行 光合 作用 的 那 部 分 太阳 
辐射 ,是 形成 植物 干 物质 的 能 量 来 源 。 作 为 一 种 重 
要 的 气候 资源 ,是 光合 潜力 潜在 产量 和 作物 生长 模 
拟 研 究 .土壤 碳 固定 模拟 研究 中 不 可 缺少 的 基础 数 
HE? 。 另 一 方面 ,光合 有 效 辐 射 还 是 水 循环 过 程 中 
生物 蒸腾 作用 的 主要 控制 因子 ,特别 是 在 生物 量 和 
作物 产量 模式 模拟 研究 中 ,光合 有 效 辐 射 更 是 不 可 
缺少 的 参数 之 一 2 。PAR 在 光合 作用 中 也 扮演 了 
重要 的 角色 中, 在 人 射 到 地 表 的 过 程 中 ,被 植被 吸 
收 和 贮存 在 生物 质 中 。 为 了 提高 对 光合 作用 和 碳 循 
环 的 模拟 ,PAR 的 估算 准确 性 非常 重要 ,其 精确 测 
量 有 助 于 提高 净 第 一 性 生产 力 (NPP ) 模 型 的 准确 
率 ,促进 C0, 源 / 汇 研究 “-”。 因 此 ,光合 有 效 辐 射 
研究 在 生态 学 、 农 学 以 及 气候 研究 中 都 具有 重要 的 
意义 和 价值 。 

尽管 光合 有 效 辐射 在 生态 系统 模型 应 用 中 极 具 


模型 的 参数 依赖 于 独立 的 辐射 站 ,而 不 同 的 辐射 站 
之 间 却 不 同 ,可 能 没有 校准 参数 的 区 域 PAR 存在 较 
KRA. 而且, 地 面 太阳 辐射 观测 站 相 较 于 其 他 
常规 的 气象 站 点 非常 稀少 ,以 地 面 太阳 辐射 站 的 
SSR 与 经 验 模型 估算 PAR 并 进行 空间 地 理 插值 获 
Hux bk KEE PAR, WH KT RAR”. A 
前 ,全 球 PAR 网 格 化 产品 空间 分 辩 率 都 比较 大 ,如 
TOMS 卫星 的 PAR p= Gh? (CERES 和 ISCCP_L 的 
PAR 产品 "等 ,空间 分 辨 率 都 超过 100 km, 而 
MODIS 只 有 海洋 PAR 产品 。 

随 着 卫星 观测 技术 的 发 展 ,以 辐射 传输 模型 为 
理论 基础 ,陆续 有 各 种 适应 特定 传感器 的 PAR 估算 
方法 提出 ,在 空间 分 辩 率 和 精度 上 取得 了 较 大 的 进 
FR) 。 在 时 间 分 辩 率 上 ,中 低 分 辩 率 的 多 光谱 遥 
感 卫 星 ,一 般 每 1 ~2 d 甚至 更 久 可 获取 一 次 全 球 综 
合 信息 ,对 于 全 球 / 区 域 的 观测 ,遥感 信息 是 一 个 瞬 
间 信 息 。 通 过 相应 的 模型 和 算法 可 估算 PAR 的 瞬 
时 值 ,而 根据 生态 模型 的 输入 要 求 , PAR 的 日 分 布 


潜力 ,但 它 并 不 是 气象 站 的 常规 观测 项 目 。 地 面 
PAR 可 由 地 面 太 阳 辐 射 (SSR，, surface solar radia- 
tion) 和 辐射 站 的 其 他 参数 估算 得 到 。 由 于 SSR 
和 PAR 存在 较 好 的 线性 关系 ,经 过 校准 的 经 验 模型 
可 以 准确 地 来 估算 该 区 域 的 PAR。 但 是 这 些 经 验 


值 更 具有 应 用 价值 "5 。 由 PAR 瞬时 值 估算 日 均值 
的 方法 较 少 ,多 数 参 考 SSR 日 均值 的 估算 方法 , 主 
要 有 以 下 几 种 :@ 加 入 间隔 时 间 , 简 单 求 均值 '”; 
D 参考 太阳 辐射 变化 规律 ,将 日 均值 与 瞬时 值 建立 
= fa BMA ~, 该 方法 在 晴朗 天 气 下 应 用 较 
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好 ,在 多 云天 气 下 效果 较 差 ;@ 将 PAR 简化 为 仅 与 
太阳 平均 天 顶 角 相关 ,该 方法 假设 当地 时 间 10:30 
代表 了 上 午 的 天 气 情况 ,13:30 代表 了 下 午 的 天 气 
情况 ;@ 根据 观测 站 多 年 数据 ,建立 瞬时 值 与 
日 均值 的 经 验 模型 2 ,该 方法 估算 的 辐射 是 SSR, 
简单 有 效 , 也 是 笔者 参考 和 改进 的 模型 与 方法 ,并 应 
用 于 PAR 日 均值 估算 。 

本 文 以 美国 国家 海洋 和 大 气管 理 局 地 面 辐 射 网 
络 的 7 个 PAR 辐射 数据 观测 站 为 研究 对 象 ,利用 
2006 一 2016 年 的 PAR 观测 数据 ,从 时 间 尺 度 和 空间 

盖 范 围 上 ,探究 PAR 瞬时 值 与 日 均值 ,以 及 PAR 
理论 值 之 间 的 估算 模型 。 为 从 瞬时 值 估算 日 均值 提 
供 一 种 简单 快捷 较为 精确 的 方法 。 瞬 时 时 间 参 考 拱 
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4% MODIS 传感器 的 Terra 和 Aqua 卫星 过 境 时 间 当 
地 时 间 一 一 10:30 和 13:30。 


1 研究 数据 和 方法 


1.1 研究 数据 来 源 

数据 来 源 于 美国 国家 海洋 和 大 气管 理 局 地 面 辐 
Bt pz (SURFRAD, Surface Radiation Budget Net- 
work) ,该 辐射 网 络 一 共有 7 个 观测 站 , 2008 年 之 前 
数据 每 间隔 3 min 记录 1 次 ,之 后 数据 每 隔 1 min 记 
录 1 次 ,单位 为 W. m“”。 该 PAR 辐射 网 络 分 布 于 
美国 7 个 州 ,主要 位 于 中 高 纬度 地 区 ,详细 的 经 纬度 
和 海拔 信息 如 表 1 所 示 。 


表 1 SURFRAD 观测 站 信息 
Tab.1 Basic information of the seven SURFRAD stations 


数据 站 纬度 (°N) ARE CE) 海拔 /m 美国 ( 州 ) 开始 时 间 
Bondville 40.06 88.37 213 伊利 诺 伊 州 1995 -01 

Fort Peek 48.31 105.10 634 蒙 大 拿 州 1995 -01 
Goodwin Greek 34.25 89.87 98 密西西比 州 1995 -01 
Table Mountain 40. 125 105.237 1 689 BLY ti Zn 1995 - 06 
Desert Rock 36. 626 116.018 1 007 内 华 达州 1998 -03 
Penn State 40.72 77.93 376 BEM JE WIN 1998 -06 
Sioux Falls 43.73 96. 62 473 南达科他 省 2003 -06 

对 地 观测 领域 中 , MODIS 数据 以 其 全 球 覆 盖 ， PAR tary = 40，SSRueos (2) 


较 高 的 时 间 分 辨 率 .多 空间 分 辨 率 .多 光谱 分 辩 率 以 
及 免费 的 信息 接受 和 产品 服务 政策 ,在 目前 研究 全 
球 尺度 下 大 气 海洋 和 陆地 生化 过 程 应 用 非常 广 
泛 所 ] 。 多 种 以 辐射 传输 模型 为 理论 基础 的 PAR ih 
算 方 法 ,数据 源 多 是 MODIS 数据 ,因此 ,本 文选 择 的 
瞬时 时 间 为 搭载 了 MODIS 传感器 的 两 颗 了 卫星 Terra 
和 Aqua 的 过 境 时 间 一 一 当地 时 间 10:30 和 13:30, 
而 实际 上 卫星 数据 获取 时 间 10:00 一 11:00 和 
13:00 一 14 :00, 因 此 ,本 文 假设 卫星 过 境 时 间 的 瞬时 
值 代 表 了 该 小 时 的 平均 情况 。 
1.2 研究 方法 与 评价 指标 

通常 情况 ,一 天 中 85% 的 日 辐射 量 在 当地 正午 
的 前 2 h 至 后 2 h SOY ,因此 ,笔者 提出 一 个 
简单 有 效 的 线性 回归 模型 ,根据 当地 时 间 10 :30 和 
13 :30 左右 的 瞬时 值 以 及 PAR 大 气 上 界 理论 值 作 
为 修正 ,来 估算 日 均值 。 


PAR gait = A, “ PAR, ning 


A, . PAR ieory +A, (1 ) 


INEP : PAR aai SPAR moming \PARwenom 分 别 为 观测 站 的 
PAR 日 均值 ,上 午 瞬 时 值 和 下 午 瞬时 值 ;PAR,,, 为 
大 气 上 界 PAR 值 ,由 大 气 上 界 太 阳 辐 射 值 SSR neo 
RE A ~ 4 为 拟 合 系数 。 

PAR 作为 太阳 辐射 的 一 部 分 , 约 占 太阳 总 辐射 
能 量 的 50% 左右 ,在 很 多 研究 中 由 地 面 接收 到 的 太 
阳 辐 射 值 和 经 验 系 数 0.48 或 0.5 直接 获得 ,不 同 的 
经 验 模型 系数 -在 0.4 ~0.5 之 间 变 动 ,在 太阳 
辐射 进入 大 气 后 , 受 云 水 汽 . 气 溶胶 等 因子 的 影响 ， 
太阳 辐射 经 吸收 .散射 后 减少 ; EEE, PAR 值 也 产生 
变化 。 考 虑 到 天 气 因 素 的 影响 ,引入 大 气 上 界 太 阳 
日 辐射 量 对 线性 回归 模型 进行 修正 , 视 为 理论 值 。 


SR = ,| Ei (3) 
SSR iheory = 2 ( T/T) < | ( sindsing + 
T J 
cos6cospcosT ) dr (4) 
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FOR 区 


PE, ERA A 1 367 Wm 一 ;ro 人 是 日 
地 相对 距离 ,与 经 纬度 相关 ;6、p Al 分 别 是 当地 的 
纬度 、 太 阳 赤 纬 和 水 平面 日 出 时 角 , 太 阳 赤 纬 和 日 出 
时 角 及 经 纬度 相关 。 

相 较 于 仅仅 用 瞬时 值 与 日 均值 建立 线性 回归 模 
型 ,引入 不 受 天 气 情况 影响 的 大 气 上 界 PAR 理论 值 
作为 修正 ,可 进一步 提高 模型 的 准确 性 ,这 也 是 笔者 
对 上 文 所 述 线 性 回归 估算 模型 的 改进 。 模 型 采用 的 
评价 指标 是 均 方 误差 RMSE .偏差 MBE 和 判定 系数 
RR ,公式 如 下 : 


N 
RMSE = J4 > (observed, — predicted, )? (5) 


N 
MBE = x (observed, — predicted, ) (6) 


R =1-> (7) 


(observed, — predicted, ) 


N - 
È (observed, - observed ) g 

N 

py 
t=1 
xP : observed, „predicted, Œ t 时刻 的 观测 值 和 模拟 
值 ;N 是 样本 个 数 。 


2 结果 与 分 析 


2.1 PAR 瞬时 值 与 日 均值 变化 分 析 

实际 上 卫星 数据 获取 时 间 在 当地 时 间 10:00 一 
11:00 和 13:00 一 14:00 ,所 以 ,以 观测 站 10:00 和 
13:00 1 min 频率 下 的 小 时 均值 代替 瞬时 值 。PAR 
观测 数据 的 日 均值 由 日 出 .日 落 时 间 内 的 1 min 频 
率 下 的 观测 值 求 得 。 

PAR 受 天 气 情 况 影 响 较 大 ,无 论 是 瞬时 值 还 是 


一 一 当地 时 间 10:30 


瞬时 值 (W m 


ChinaXiv 合 作 期 刊 
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日 均值 ,日 与 日 之 间 差 异 均 较 大 ,但 总 体 呈 现 出 冬 春 
季 低 、 秋 夏季 高 的 特征 ,季节 变化 明显 。 如 图 1 所 
示 , 以 Bondville 观测 站 2016 年 PAR 的 瞬时 值 与 日 
均值 为 例 ,PAR 日 均值 与 瞬时 值 均 波动 较 大 ,上 午 
和 下 午 的 瞬时 值 总 体 处 于 50 ~ 400 W.m ,日 均 
值 总 体 处 于 20 ~ 160 W- m”, RF BEHE, 
此 ,以 瞬时 值 的 天 气 情况 代替 上 午 、 下午 的 天 气 情 
况 ,可 能 带 来 一 定 误差 。 全 年 PAR 变化 呈现 出 正弦 
函数 的 变化 规律 ,其 他 观测 站 日 均值 年 变化 类 似 。 
2.2 观测 站 数据 模拟 结果 验证 与 分 析 

光合 有 效 辐射 主要 受 地 理 位 置 海拔 和 天 气 情 
况 影响 ,基于 辐射 传播 模型 的 遥感 估算 方法 考虑 了 
上 述 影 响 因素 估 算出 瞬时 值 。 在 此 基础 上 ,笔者 探 
究 由 瞬时 值 估算 日 均值 的 线性 回归 模型 是 否 受 地 理 
位 置 海拔 以 及 气候 因素 的 影响 。 同 上 文 所 述 ,以 观 
测 站 10:00 和 13 :00 在 1 min 频率 下 的 小 时 均值 代 
替 瞬 时 值 。PAR 观测 数据 日 均值 由 日 出 .日 落 时 间 
内 的 1 min 频率 下 的 观测 值 求 得 。 

对 每 年 的 数据 进行 拟 合 后 发 现 ,同一 个 观测 站 
不 同年 份 的 拟 合 系数 非常 接近 , 故 认为 该 线性 回归 
方程 尽 可 能 减少 了 气候 影响 ,因此 ,对 同一 个 数据 站 
的 拟 合 系数 采用 多 年 拟 合 系数 的 平均 值 。 利 用 11 a 
观测 数据 对 拟 合 结果 进行 验证 ,图 2 分 别 是 7 个 站 
点 的 模型 拟 合 结果 ,RMSE 和 MBE fy AA et W 

从 图 2 可 以 看 出 ,此 线性 回归 方程 应 用 于 不 同 
区 域 的 单个 光合 有 效 辐 射 观测 站 ,RMSE MBE il R? 
均 较 小 ,其 中 RMSE 均 低 于 10 W-m’?, 均 大 于 
0.96 , 拟 合 结果 较 好 , 拟 合 系数 如 表 2 所 示 。 


当地 时 间 13:30 
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日 均值 (W -m 
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图 1 Bondville 观测 站 2016 年 PAR 瞬时 值 与 日 均值 变化 


Fig.1 Changes of PAR instantaneous and daily values at Bondville Station in 2016 
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(b)Fort Peck 
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MBE=0.324 7 
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(f)Penn State 
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图 2 2006 一 2016 年 7 个 观测 站 模拟 
Fig.2 Validation of the simulated results at seven SURFRAD stations during the period from 2006 to 2016 
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R2 数据 站 拟 合 系数 
Tab.2 Fitting coefficients for the SURFRAD stations 


数据 站 拟 合 系数 
A, A, Ag XA; A, 
Bondville 0.163 727 273 0. 146 281 818 0.077 772 727 — 19.573 545 5 
Fort Peek 0. 162 254 545 0. 149 872 727 0.085 945 455 -— 18.340 763 6 
Goodwin Greek 0.163 7 0. 141 018 182 0.095 127 273 — 22.742 227 3 
Table Mountain 0. 164 154 545 0. 126 881 818 0.0909 454 55 — 19.496 663 6 
Desert Rock 0.150 118 182 0.141 8 0. 122 045 455 — 24.245 290 9 
Penn State 0. 157 990 909 0. 149 772 727 0.069 427 273 -17.026 554 5 
Sioux Falls 0. 1573 090 91 0.160 418 182 0.078 827 273 - 19. 125 672 7 
平均 0. 159 893 506 0.145 149 351 0.088 584 416 — 20.078 674 
#3 SURFRAD 观测 站 模拟 结果 
Tab.3 Simulated results at the seven SURFRAD stations 
BON FPK GWN TBL DRA PSU SXF 全 部 
RMSE/(W m~?) 7.707 9 8.696 7 8.310 2 9.001 8 7.218 5 8.516 1 7.742 6 8.863 1 
MBE/(W - m`?) -0.030 01 0.324 65 -3.128 1 0.074 703 -0.767 04 0.0193 02 —0. 008 38 0.433 57 
R? 0. 980 6 0.978 5 0.980 6 0.968 6 0.983 8 0.975 9 0.981 1 0.976 97 


R°=0.977 
RMSE=8.863 1 
MBE=0.433 6 
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图 3 7 个 观测 站 2017 年 模拟 结果 验证 
Fig.3 Validation of simulated results at seven SURFRAD 


stations in 2017 


从 表 2 可 以 看 出 ,位 于 不 同 地 理 位 置 的 多 个 观 
测 站 拟 合 系数 也 较 接近 ,存在 微小 的 差异 。 为 探究 
该 回归 模型 统一 后 的 拟 合 系数 是 否 适 用 于 不 同 区 
域 ,笔者 利用 全 部 数据 站 多 年 拟 合 系数 的 平均 值 作 
为 最 终 的 拟 合 系数 ,应 用 于 该 模型 中 ,输入 2017 年 
7 个 观测 站 上 午 、 下 午 瞬 时 值 数据 估算 日 均值 ,并 用 
观测 值 进行 验证 ,多 个 数据 站 模拟 验证 结果 如 图 3 
所 示 ,误差 分 析 如 表 3 所 示 , RMSE 和 MBE 的 单位 


-2 
是 W'm 。 


从 表 3 可 以 看 出 ,该 模型 和 固定 的 拟 合 系数 应 
用 于 多 种 经 纬度 ,海拔 高 度 和 气候 条 件 下 的 观测 站 , 
拟 合 效果 仍然 较 好 ,误差 较 小 , RMSE 为 8. 863 1 
W +m ,及 为 0.976 97, 根据 卫星 遥感 数据 估算 
的 瞬时 PAR 值 ,可 以 考虑 应 用 此 模型 和 系数 ,由 上 
午 \ 下 午 瞬 时 值 估 算 日 均值 。 
2.3 观测 站 数据 模拟 结果 验证 与 分 析 

与 其 他 方法 相 比 ,利用 间隔 时 间 求 均值 所 方 
法 ,其 在 晴朗 天 气 下 的 RMSE 为 32. 621 03 W - 
m? ZARA FHI RMSE 为 63.785 53 W m”, £ 
种 天 气 条 件 下 的 RMSE X 45.7193 W- m°; iit 
建立 三 角 函 数 关系 由 瞬时 值 估算 日 均值 的 方法 "”， 
其 在 多 个 区 域 RMSE 范围 为 33.6~147.1 W+ m°; 
仅 利 用 上 午 和 下 午 瞬 时 值 建立 估算 模型 的 方法 7， 
JE RMSE 为 53.33 W nm- 一 ;前 3 种 方法 的 瞬时 值 由 
卫星 数据 估算 而 来 ,因为 估算 的 瞬时 值 本 身 存在 一 
定 误差 ,所 以 ,由 这 样 的 估算 值 作为 瞬时 值 输入 ,在 
估算 日 均值 时 比 利 用 观测 站 的 瞬时 值 作为 输入 误差 
可 能 会 偏 大 。 

为 了 验证 笔者 所 提出 的 线性 回归 模型 以 基于 卫 
星 数据 估算 的 瞬时 值 为 输入 ,估算 日 均值 是 否 准 确 
有 效 ,基于 辐射 传播 模型 "和 MODIS 卫星 数据 估 
算 了 上 午 和 下 午 的 瞬时 值 。 该 辐射 传播 模型 由 唐 文 
REE 2017 年 提出 ,考虑 多 种 气候 条 件 因素 ,整体 
误差 较 小 。 验 证 数据 为 观测 站 过 境 时 间 小 时 平均 数 
据 ,验证 结果 如 图 4 所 示 , RMSE 和 MBE 的 单位 是 
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图 4 Bondvile 观测 站 范围 基于 MODIS 数据 估算 的 
瞬时 值 验证 结果 
Fig.4 Validated results of the instantaneous values estimated 
by MODIS data at Bondville Station during the period 
from 2015 and 2016 
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图 5 基于 估算 的 瞬时 值 和 线性 回归 模型 估算 的 

日 均值 验证 结果 


Fig.5 Validated results of the estimated instantaneous values 


and the daily values based on the linear regression model 


Wem”, 

如 图 4 所 示 , 由 卫星 数据 和 辐射 传播 模型 估算 
的 上 午 、 下 午 瞬 时 值 存在 一 定 的 估算 误差 ,但 误差 尚 
在 合理 范围 内 。 因 为 PAR 的 影响 机 制 复杂 ,一 方面 
误差 来 源 于 MODIS 数据 估算 PAR 时 所 采用 的 方法 
本 身 存 在 的 误差 ; 另 一 方面 来 源 于 MODIS 多 个 数据 
产品 ,在 云天 气 情况 下 多 种 产品 存在 大 量 缺 失 值 , 需 
要 通过 其 他 数据 如 MODIS 的 月 均值 和 美国 国家 环 
境 预 报 中 心 (NCEP) 再 分 析 资 料 进行 插值 填充 。 估 
算 的 卫星 瞬时 值 作 为 本 文 提 出 的 线性 回归 模型 的 输 
人 来 估算 日 均值 ,探究 在 输入 的 瞬时 值 存在 一 定 误 
差 的 情况 下 ,结合 上 午 、 下 午 瞬 时 值 ,引入 理论 值 进 
行 修正 后 ,验证 本 文 提出 的 线性 回归 模型 效果 。 验 
证 数据 为 观测 站 日 均值 ,模拟 结果 如 图 5 所 示 ， 
RMSE 和 MBE HAIE W m°’ 

从 图 4 和 图 5 可 以 看 出 ,即使 在 估算 的 瞬时 值 
存在 一 定 误 差 的 情况 下 ,利用 本 文 提出 的 线性 回归 
模型 和 拟 合 系数 依旧 很 好 地 估算 出 了 较为 准确 的 日 
均值 ,RMSE X 14.274 6 W - m°, R? 为 0.926 82, 
而 上 文 所 述 的 将 PAR 简化 为 仅 与 太阳 平均 天 项 角 
相关 求 日 均值 的 方法 " ,其 输入 的 瞬时 值 由 卫星 数 
据 和 辐射 传播 模型 估算 而 来 时 ,RMSE 为 14.6 W- 
m ,MBE 为 3.7 W- m”, R 为 0.92, 与 本 文 模型 
的 验证 结果 基本 持平 ,但 考虑 到 其 应 用 的 是 2008 年 
与 2009 年 数据 ,而 且 在 辐射 传播 模型 中 及 其 所 用 遥 
感 数 据 本 身 存在 误差 ,所 以 ,误差 总 体 来 源 于 卫星 估 
算 的 瞬时 PAR; 其 输入 的 瞬时 值 来 源 于 观测 站 时 ， 
RMSE 为 11.9W-.m ,MBE 为 2.7W.m ,R 为 
0.94 ,总 体 低 于 本 文采 用 的 线性 回归 模型 估算 方法 
结果 ,RMSE 为 8.863 1 W - m”, MBE 为 0.433 57 
W m°, R 为 0.976 97。 由 此 可 见 , 本 文 提出 的 线 
性 回归 模型 是 可 靠 的 ,可 以 简单 ` 有 效 .准确 地 根据 
瞬时 值 估算 PAR 日 均值 。 


3 结论 


本 文 基于 SURFRAD 地 面 辐 射 观测 网 络 中 不 同 
经 纬度 分 布 海拔 高 度 下 的 7 个 观测 站 2006 一 2016 
年 PAR 数据 ,分 析 了 一 年 内 的 PAR 日 均值 变化 和 
多 年 年 均值 变化 ,日 均值 在 20 ~160 W:m”? BR 
位 于 高 海拔 的 TBL 观测 站 年 均值 在 100 ~ 105 W - 
m~’ ,其 他 观测 站 年 均值 在 62 ~85 W + m”, 波 动 较 
小 。 可 以 看 出 PAR 日 均值 变化 较 大 ,年 均值 变化 较 


FOR 区 


平稳 ,但 是 季节 规律 明显 , 冬 春 季 低 、 夏 秋季 高 。 

以 辐射 传输 模型 为 理论 基础 ,利用 遥感 卫星 数 
据 估算 瞬时 PAR 的 方法 逐渐 成 熟 。 由 瞬时 值 估算 
日 分 布 值 , 笔 者 提出 一 种 简单 ` 有 效 的 线性 回归 模 
型 ,根据 PAR 瞬时 值 以 及 大 气 上 界 SSR 值 估算 PAR 
日 均值 ,效果 良好 。 单 个 观测 站 的 模拟 结果 RMSE 
值 均 低 于 10 W m? RANS W m eR? 均 
高 于 0. 96 ,最低 为 0.968 6。 应 用 该 方法 对 2017 年 
7 个 观测 站 数据 进行 估算 并 验证 ,模拟 结果 显示 
RMSE 为 8.863 1 W .m- ,MBE 为 0.433 57 W ， 
m`’, R? 为 0.976 97, 

本 文 建立 的 估算 模型 仅 根 据 PAR 上 午 、 下 午 的 
瞬时 值 及 与 跟 经 纬度 相关 的 大 气 上 界 SSR 理论 值 
进行 估算 ,简化 了 模型 。 模 拟 结果 显示 ,该 模型 具有 
良好 的 可 行 性 和 可 靠 性 ,在 由 瞬时 值 估算 日 均值 的 
过 程 中 不 考虑 地 形 .经 纬度 分 布 和 和 气候 条 件 的 影响 ， 
可 以 考虑 大 范围 应 用 。 

综 上 所 述 ,本 文 基于 SURFRAD 地 面 辐射 观测 
网 络 数据 建立 了 线性 回归 模型 ,由 PARR 上午、 下 午 
瞬时 值 估算 日 均值 ,避免 再 次 考虑 地 理 区 域 . 气 候 情 
况 和 海拔 高 度 等 因素 的 影响 ,减少 了 传统 方法 的 误 
差 。 无 论 是 基于 观测 站 的 瞬时 值 数据 ,或 者 基于 辐 
射 传播 模型 和 卫星 数据 估算 的 瞬时 值 数据 , 均 可 以 
准确 有 效 地 估算 出 日 均值 。 为 后 续 大 范围 基于 遥感 
数据 估算 瞬时 PAR 值 后 估算 日 均值 提供 了 可 行 可 
靠 的 方法 。 结 合 遥 感 数据 估算 PAR 日 均值 对 研究 
植被 碳 循环 具有 重要 的 意义 ,有 助 于 进一步 提高 植 
被 净 生 产 力 估算 模型 的 准确 性 。 
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Estimation of Daily Value from PAR Instantaneous Value 


ZHANG Jing'*, XIE Ya-nan', WANG Ming-quan’, WANG Mao-hua” 
(1. Key laboratory of Specialty Fiber Optics and Optical Access Networks , Shanghai Institute of Advanced Communication 
and Data Science ,Shanghai University ,Shanghai 200000 , China ; 
2. Shanghai Carbon Data Research Center ,Shanghai Advanced Research Institute , Chinese Academy of Sciences , 
Shanghai 200000 , China) 


Abstract: — Photosynthetically active radiation (PAR) commonly refers to solar radiation in the visible (400 — 700 
nm) part of the electromagnetic spectrum. Despite its potential and critical application in terrestrial ecosystem mod- 
elling ,surface PAR is not typically a routine observation at meteorological station. Remote sensing can be used to 
derive the spatiotemporally continuous instantaneous PAR, but the researches about daily PAR estimation are defi- 
cient,and their accuracy cannot be guaranteed. In this paper,a linear regression model was developed to calculate 
the daily PAR based on Surface Radiation Budget Network (SURFRAD) data through the instantaneous and theo- 
retical PAR. The results indicated that the model was effective with root mean square errors (RMSE) about 8.863 1 
W +m’ with coefficient of determination (R*) about 0.977 without the limitations of areas. 
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